Non-holonomic Quantum Devices 
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We analyze the possibility and efficiency of non-holonomic control over quantum devices with 
exponentially large number of Hilbert space dimensions. We show that completely controllable 
devices of this type can be assembled from elementary units of arbitrary physical nature, and can be 
employed efficiently for universal quantum computations and simulation of quantum field dynamics. 
As an example we describe a toy device that can perform Toffoli-gate transformations and discrete 
Fourier transform on 9 qubits. 
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I. INTRODUCTION 



What is the difference between a classical and a quan- 
tum device? Clearly it is not in the physical laws gov- 
erning their dynamics, since Classical Mechanics follows 
from Quantum Mechanics as a limiting case, when me- 
chanical action for each degree of freedom is much larger 
than the Planck constant fi. Hence, all classical devices 
are quantum as well, and the basic difference between 
them is rather in the quantities of interest and in the 
interactions under control. Typically, the operators of 
main physical quantities have smooth dependence of their 
semiclassical matrix elements on the indices numerating 
the energy eigenstates, and therefore a state of the device 
is characterized by the position of the center of the Ehren- 
fest wave packet in phase space. The average quanti- 
ties are determined as functions of this position, whereas 
the finite packet width results in uncontrolled "quantum 
noise" and is considered as an obstacle for the correct 
operation of the classical device in the quantum limit. 

The situation is different in the essentially quantum 
limit, where the action for each degree of freedom is of 
the order of h. Then, the matrix elements are not smooth 
anymore, and the consistent description of an iV-level de- 
vice relies not only on quantum averages of operators, but 
also on all their higher moments as well. Such descrip- 
tion requires exhaustive information about the state of 
the system, as given by a vector in the iV-dimensional 
Hilbert space of the system. Building a completely con- 
trolled quantum device in practice implies control over 
all the moments and therefore is a challenging task. It 
promises, however, adequately important practical ben- 
efits: coherent control of molecules, quantum cryptogra- 
phy, and quantum computation are some of the potential 
applications fl|-|l5[|. 

In this paper we describe a scheme for constructing 
completely controllable quantum devices. We show that 
quantum systems perturbed in a certain time-dependent 
way become "non-holonomic" , which means that as a 
result of the perturbation all global constraints on the 



dynamics are removed and the system becomes fully con- 
trolled (Sec. ||). We then describe a simple, completely 
controllable "unit cell" that can serve as a building block 
for compound devices of arbitrary size (Sec. Ill), and 
show in particular that it can implement the Toffoli gate 
(Appendix [X]) . We give examples of compound devices 
that can be employed efficiently for universal quantum 
computations and simulation of quantum field dynam- 
ics (Sec. IV). Finally, we describe a toy device that can 
perform quantum computations on 9 qubits and show in 
particular how it can perform the discrete Fourier trans- 
form on 9 qubits (Sec. [v]). 



II. NON-HOLONOMIC CONTROL 

The idea of controlling a system by forcing it to 
have globally unconstrained — non-holonomic — dynamics 
is natural, since in order to ensure an arbitrary evolution 
one has first to get rid of the restrictions posed by the 
existing integrals of motion and all other constraints. In 
the non-holonomic control scheme, the system evolution 
is determined by an unperturbed Hamiltonian Hq and 
a number of perturbations C,-Pj of fixed operator struc- 
ture Pi and controllable strengths Ci that are applied to 
the system, so that the evolution is given by the time- 
dependent Hamiltonian 



H{t) = H {) + Y,C l (t)Pi. 



(1) 



The system becomes non-holonomic and completely con- 
trollable if the commutators of all orders of Ho and the 
Pi span the space of Hermitian operators in the Hilbert 
space of the system, that is, if an arbitrary Hermitian 
operator can be represented as a linear combination of 
the operators 



H , Pi, [H , A], [PuPj], [P,, [Pj,P k 



(2) 



Note that at most N 2 linearly independent terms are 
needed for an N dimensional Hilbert space. 
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The control scheme consists of two steps: (i) verifi- 
cation that the perturbations induce non-holonomic dy- 
namics; and (ii) finding particular time dependencies for 
the perturbations that effect a given desired control. Step 
(i) is straightforward — by inspecting the commutation re- 
lations between the explicitly written Hamiltonian and 
perturbation operators, one checks if the system under 
consideration is indeed non-holonomic. But step (ii) re- 
quires more art — one has to put the system in such condi- 
tions that all unwanted outcomes, present in abundance 
in a system with no constraints, experience a destructive 
quantum interference. 

Let us consider a quantum system of N = 2™ levels, 
composed of n interacting two-level subsystems. To be 
specific we speak about two-level atoms in a laser field, 
although it could as well be any other quantum object, 
such as interacting spins in a magnetic field, Josephson 
junctions, Rydberg atoms, rotating molecules, quantum 
dots on a surface, etc. The only requirement is that 
the object must be subjected to a non-holonomic con- 
trol, since only in this case it can perform any desired 
operation, no matter what the physical interactions in 
the system are. The choice of a practical realization of 
a non-holonomic system will therefore depend mainly on 
optimization of technical parameters such as simplicity 
and cost-effectiveness of construction, lifetime of quan- 
tum coherence ([3], precision of available controls, and 
so on. 

A crucial issue that determines the strategy of con- 
struction is the required extent of immediate universality 
of the control. In principle, one can think about complete 
and direct physical control over a 2 n -level quantum sys- 
tem, even for a large n, which implies the ability to ensure 
an arbitrary evolution of the system, given by any prede- 
termined 2" x 2" unitary matrix U, and which requires 4™ 
physical control parameters. For this purpose one should 
find an algorithm that determines these controls for any 
given U. ft might be difficult to find such algorithm, 
and even if found, its application will require an enor- 
mous computational work that grows exponentially with 
n, and will therefore be intractable. In addition, the cost 
of physically implementing the huge number of 4™ control 
parameters seems too high a price to pay for this kind 
of universality, which may not even be needed for practi- 
cal purposes. For these reasons, one should presumably 
give up direct universality and search for specialized ways 
to build quantum devices for each particular task, with 
number of controls that is not exponentially larger then 
what is specifically needed. 

III. COMPLETELY CONTROLLED UNIT CELL 

One way to construct a completely controlled but not 
immediately universal quantum device is to build it up 
from small parts, "unit cells" , each of which is non- 
holonomic and therefore directly and universally control- 
lable. The proper functioning of the device relies then on 
the appropriate connection of the cells p~7[ . In this way 



the universality of the device is obtained indirectly, not 
by applying a huge number of controls, but by smartly 
connecting the cells and choosing the sequence of oper- 
ations performed. There is no general prescription how 
to construct a particular device; this requires expertise 
in the art of "programming" the operations of the cells 
and their interactions. 

A. Cell structure 

An example of a completely controlled unit cell is 
shown in Fig. [l| It consists of three two-level atoms, 
each with ground and excited states |0) and |1), hav- 
ing distinct transition frequencies w", , and lu 3 . The 
atoms have dipole-dipole interaction between themselves 
and are coupled to two external fields: an electromag- 
netic field E u — £ u cos Lot of nearly resonant frequency 
u), and a static electric field Es- The dipole-dipole inter- 
action is fixed and determines the principal, unperturbed 
Hamiltonian of the system, Hq, while the external fields 
provide two controllable perturbations, P u and Ps- The 
Hilbert space of the system has a "computational basis" 
of N — 2 3 = 8 states, \x) = \x 2 xix ) = \x2)\xi)\xq) , 
x = 0, 1, . . . , 7, where the state of the ith atom encodes 
the ith binary digit of x = ^2 r —oX r 2 r as a qubit [see 
Fig. [j](b)]. The crucial requirement is the non-holonomic 
character of the interaction. It implies that Hq, P u , Ps, 
and their commutators of all orders span the linear space 
of 8 x 8 Hermitian matrices |jgj] . This is indeed the 
case for the system shown in Fig. [j], which has princi- 
pal Hamiltonian and perturbations given, in the compu- 
tational basis and assuming resonant approximation, by 
the matrices 
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Here = didj/R^ is the dipole-dipole coupling of the 
ith and jth atoms at distance Rij , with di the ith atom 
dipole matrix element, and Vi — E^di is the dipole cou- 
pling of the ith atom to the external electromagnetic 
field. The excitation energy detunings of single atoms 
Ai = h(ujf — oj) determine the detunings of pairs of atoms 
Aij — Ai+Aj and the total detuning A a = + 
Their values can be changed by variation of a static elec- 
tric field Es (Stark effect), which results in energy shifts 
Ai — > Ai + Ai for single atoms, where A, = cuEs depend 
on atom-specific electric permeability constants Qfj, and 



similar shifts A,-, 



A,- and A a = Ai + A 2 + A 3 for 



two and three atomic detunings respectively. 

Note that by a proper choice of Aj and u> one can set 
two of the three Ai to zero. Moreover, to simplify the pre- 
sentation we also set to zero the third Ai, which would 
otherwise remain just a part of Hq. Hence, hereafter all 
Aj denote just the deviations from zero resulting from 
the variation of the Stark field Es- The latter together 
with the amplitude £ w serve as time dependent control 
parameters, C's and C w respectively. The matrices P$ 
and P u contain therefore only the permeabilities and 
the dipole moments di respectively. 



B. Cell control 

To exert direct universal control over the unit cell we 
proceed as follows, (i) We fix N 2 — 64 consecutive time 
intervals of equal duration T in which the two pertur- 
bations will be applied to the system in an alternating 
sequence: in the fcth interval the perturbation is Pk = Ps 
for odd k and Pk — P u for even k, where k = 1, 2, ... , 64. 
The strength of Pk is a controllable parameter, which we 
take to have a constant value Ck, and which denotes ei- 
ther or Eg, during the fc-th time interval, depending 
on the parity of k Jy^. Thus, the system evolution is 
given by a Hamiltonian which is constant in each inter- 
val: 



H(t) = H + C k P k te[(k-l)T,kT]. 



(6) 



(ii) We find 64 positive C k values for which the total evo- 
lution of the system will be the identity transformation: 



64 



U(t = 64T) = J| exp 



fc=i 



-(H + C k Pk)T 



(7) 



To this end, we first solve the "8th root" of Eq. (^) , 
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FIG. 1. Realization of a unit cell: A compound system of 
three two-level atoms interacting with external electromag- 
netic and static electric fields, (a) The ith atom has ground 
and excited states \0)i and \ with excitation energy Ai + Ai 
that can be modified by the static field; transition amplitude 
in the electromagnetic field is Vi] the dipole-dipole coupling 
of the ith and jth atoms is Dij. (b) The computational ba- 
sis states and their relation to matrix elements of the prin- 
cipal Hamiltonian Hq and the perturbations and Ps of 
Eqs. (||). 



U(t = 8T) ee Y[ exp 



fc=i 



■(H + C k Pk)T 



JV8, (8) 



by minimizing the coefficients of the characteristic poly- 
nomial of U(t = 8T) fl4| . This gives a sequence of posi- 
tive values, Ci, C2, . . . , Cs, for which U(t — 8T) has the 
eigenvalues e 27rmi / 8 , m — 1, 2, . . . , 8, and hence satisfies 
[U(t = 8T)] 8 = I nondegenerately. Repeating this se- 
quence 8 times we obtain the required 64 Ck- (iii) Now, 
by small variations SCk of the Ck values we can obtain 
any unitary transformation t/ e in a small neighborhood 
of the identity transformation: 



1)4 



U(t = 64T) ee J| exp 



fc=i 



-(H + [C k + SCk}Pk)T 



= U f . 



(9) 



Indeed, we can present this "small" transformation as 

U e = exp{-iHe), (10) 

with dimensionless 8x8 Hermitian Hamiltonian TL which 
is bounded as < 1 and is multiplied by a small pa- 
rameter e > 0. Now the variations 5C k are determined 
to first order in e by the linear equations 



64 

E 

k=l 



dU(t = 64T) 
dCk 



SC k 



-iHe. 



(11) 



Moreover, when e is sufficiently small, iterative Newton 
method refinements of the SC k yield U(t = 64T) = U e 
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with utmost accuracy p4| . (iv) Finally, to perform an 
arbitrary unitary transformation we again present it as 
U e in Eq. ([lO]), but now the parameter e may take any 
value in [0, 2ir] and will not necessarily be small. We ef- 
fect the U e by dividing it into "small" steps: we apply the 
transformation U(t = 64T) = U e / m repeatedly m times, 
with m big (e/m small) enough to allow direct control, 
and obtain 



U(t = MmT) = (U t/m ) rn = U, 



(12) 



We note that, for the problem under consideration, direct 
control is typically attainable with m < 16. Moreover, 
we expect that m = 1 will be sufficient with more power- 
ful numerical methods for solving Eq 
contrast with earlier control schemes 
unitary transformation is effected within a few control 
cycles, with accuracy that depends in principle only on 
the physical precision of the controls. 
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FIG. 2. (a) Control parameters Ck for the identity trans- 
formation /. Variations SCk effecting on the cell the transfor- 
mation U f /s, with U c = (U e / S ) 8 equal to: (b) the permutation 
U p i2', (c) the permutation U P 23; (d) the Toffoli-gate transfor- 
mation Uroff- (e) Variations SCk effecting the conditional 
phase shift B((j>) = exp(— itftHs), at </> = 7r/32, employed in 
the quantum discrete Fourier transform. 

In Fig. H we show examples of unit cell control, where 
appropriately chosen parameters Ck and variations SCk 
effect unitary transformations on the unit cell: the 
Toffoli-gate transformation (see Appendix |a|), two-qubit 
permutations Pij\cL)i\b)j — \b)i\a)j (a, b — 0,1), and the 
conditional phase shift employed in the quantum discrete 
Fourier transform (discussed in Sec. |V|). The transforma- 
tion is achieved either directly (m = 1) or by 8 repeti- 



tions (m = 8). The operators Hq, and Ps are chosen 
with arbitrary realistic values. We take Dyi = 1.1E U , 
L>23 = 0.946£;„, Di 3 = 0.86K, and T = 250h/E u , 
where E u ~ 10~ 18 erg is the typical energy scale. For 
odd k we switch off the external electromagnetic field, 
^1:2:3 = 0, and tune the atomic excitation energies by 
the Stark field E s such that Ai ;2;3 = (0.1; 0.11; 0.312)^ u . 
For even k we set Es — 0, that is Ai-2-3 = 0, and take 
Vi ;2;3 = (0.3; 0.33; 0.24)£„- 



IV. COMPLETELY CONTROLLED QUANTUM 
DEVICES 

Once completely controlled unit cells can be con- 
structed, a compound device can be assembled from such 
elements. To be efficient, the architecture of the device 
will depend on the specific function it should perform. 
In Fig. |3| we show two possible arrangements of unit cells 
for special purpose devices: the first arrangement suites 
more the purpose of quantum computing, while the sec- 
ond is more useful for simulating lattice quantum field 
dynamics. 
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(a) 



# © # © # © # © # © # # 

(b) 

# © # © # # © # © # 0. 

FIG. 3. Two possible arrangements of cells for special 
purpose devices: (a) tree-like structure for quantum compu- 
tation; (b) planar lattice for simulating dynamics of quantum 
fields. The circled numbers denote the rank of joints of the 
tree (a) or specify the order in which atoms are grouped into 
triads (b). The arrows show state exchange to parent joints. 

The first device (Fig. |(a)) is organized in a tree-like 
structure, where the quantum state of one atom in each 
cell can be exchanged with the state of an atom at the 
closest parent joint of the tree. The simplest way to make 
the exchange is to displace the atom to the parent joint, 
however, the exchange or transport of the state without 
moving the atom can be more practical. The tree-like ar- 
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chitecture and the possibility to perform all the unitary 
transformations, including all the permutations, in each 
unit cell allow one to put together and make interfering 
the states of any three two-level atoms of the device after 
at most s = 6 log 3 n state exchanges, by moving them to- 
ward the root of the tree to a common cell. Placing the 
new states back (if needed) requires the same number of 
inverse exchanges. This is a very modest number, s ~ 40, 
even for a rather large device of n ~ 10 3 with Hilbert 
space of N = 2" ~ 10 300 dimensions. Hence, all basic 
operations of quantum computation can be performed on 
any physical system comprised of non-holonomic triads 
of two-level subsystems in a tree-like structure, and each 
operation can be completed within 64 x 16 X 12 x log 3 n 
control intervals T. Note that the unity transformation 
should be applied to all other cells to preserve their states 
during the operation. 

The second arrangement of cells (Fig. |](b)) is intended 
mainly for emulating the dynamics of quantum fields on 
lattices. Of course, it can also perform general operations 
on any triad, but for a higher cost of s = 0(n l l 2 ). In this 
arrangement, after each control period of 64T the clos- 
est neighboring atoms are differently regrouped in triads 
(cells), with the original grouping repeating itself after 
three consecutive periods. Therefore, at each moment 
the change of the cell state depends on the states of the 
neighboring cells, as it should be in order to emulate the 
dynamics of the fields. Immediate analogy to the Ising 
model emerges when we restrict ourselves to small values 
of e where terms of order e 2 are negligible, and then each 
64T period plays the role of the time increment At = e. 
The evolution of such device is determined by three sums 
of effective cell Hamiltonians, = ^2, q "Hq,p, onc f° r 

each period p = 1, 2, 3, where Ti. q , p is the effective Hamil- 
tonian of the gth cell at the pth period. 

We can cast the cell Hamiltonians to sums of tensor 
products of Pauli matrices a l a , where the Greek index 
a = x,y, z denotes the matrix type and the Latin index 
i specifies the two-level atom on which it acts. Since the 
cells are under complete control, the coefficients of this 
development can be made an arbitrary function of the 
time r, and hence the effective Hamiltonian reads 

H eff (r) = Af(T)a i a + B^ j) (r)a i a ^ 

+ Cf2 k) {T)dl^b% (13) 

with implicit summation over repeated indices, where 
and (i,j,k) indicate pairs and triads of distinct 
atoms that are periodically grouped in a common cell. 
This Hamiltonian results in the evolution equation for 
the Heisenberg operators (r) , 

n d^rl = ^. (t) ^ (t) + fl^ ffc) (r)^(r)fi*(r) 

+ CiS, (r)^(r)^(r)4(r), (14) 

where the coefficients A,B,C are determined by A, B, C 
and the commutation relations of the Pauli matrices. By 



a proper choice of the coefficients A, B, C through the 
appropriate control sequences, one can simulate different 
linear and non-linear lattice models of quantum fields 
with time dependent parameters. 



V. TOY DEVICE 




FIG. 4. A toy device, composed of 9 Rydberg atoms, that 
can perform quantum computations on 9 qubits. Each atom 
is a two-level system shown schematically by double orbits. 
Atoms of different triads are excited to distinct pairs of Ryd- 
berg states. Each triad p is controlled by an external field of 
distinct frequency u) p . One atom in each triad can be excited 
to a pair of higher Rydberg states, thus forming a higher-level 
triad: (3,6,7). These excitations (depicted by arrows) corre- 
spond to state transportations. 

We now describe a toy device that can perform quan- 
tum computations on 9 qubits. An ensemble of 9 differ- 
ent Rydberg atoms is placed in a magneto-optical trap 
at low temperature, as illustrated in Fig. [|. By differ- 
ent atoms we mean atoms of different elements or iden- 
tical atoms that are excited to distinct pairs of Rydberg 
states. The best candidates for such a device are the 
long-living states corresponding to large angular momen- 
tum. By placing all the atoms in a static electric field 
one lifts the degeneracy of the magnetic quantum num- 
ber and performs tuning if needed. All the atoms ex- 
perience the dipole-dipole interaction L>ij = didj(R~j 3 ), 
where the cube of the inverse distance between atoms is 
averaged over their translational quantum states. Note, 
however, that only for almost resonant atoms this inter- 
action is important. By a proper choice of the atomic 
states and the static field Eg, we obtain three triads, 
p = 1,2,3, each comprised of three almost resonant two- 
level atoms with transition frequencies centered around a 
distinct frequency lu p . For each triad p, the interactions 
Dij give the principal Hamiltonian, while a microwave 
field E w at the frequency oj p serves as a control per- 
turbation. Transportation of the state of one atom in 
each triad to the parent joint can be performed by dipole 
or Raman 7r transitions from the initial pair of Rydberg 
levels to a higher pair. With these higher pairs assumed 
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nearly resonant with a frequency W4, atoms 3, 6 and 7 
form a higher-level triad — the parent joint of the first 
three triads — which is controlled by a forth microwave 
field E UJ4 of frequency 0J4 . 

As an example of implementing quantum computa- 
tion in the toy device, using our non-holonomic control 
scheme, we show how to perform the discrete Fourier 
transform modulo TV = 2 9 = 512 pC|j . This is the uni- 
tary transformation on 9 qubits that is given by 



F N \x) 



1 



N-l 



^2 e*p(2irixy/N)\y), 



(15) 



v=o 



where \x) and \y) are states of the system computational 
basis. The computational basis states are defined as 



F8/9 ■ 



Xl>2 £ )l, 



(16) 



with X = J2r=0 



0, 1, . . . ,N - 1 (x r = 0, 1), where 



I )i denotes the state of the ith atom — the ith qubit. The 
algorithm we employ to perform the Fourier transform is 
based on constructing the exponent in Eq. (|l^) as 

8 r 

exp(2irixy/2 9 ) = [] II cxp(i7r<y s /2^ s ), (17) 

r=0s=0 

where x' r = x 8 _ r . We begin by reversing the order in 
which the bits of the input x are stored in our 9-qubit 
register, that is, we effect the unitary transformation 



|x 8 ) 9 . . . |xi) 2 |a;o)] 



\xo)g ■ ■ • \x7) 2 \x S )l 



(18) 



by applying a sequence of state exchanges |21j. Then 
we complete the transform in 9 steps: (i) We "split" the 
first qubit (the state of atom 1) by applying the unitary 
transformation 



.4 



V2 \l 



1 



exp[ 



— ITT 

V8 



1-V2 1 

1 -1 - 



(19) 



which maps |0) -> ^(|0) + |1» and |1) - i(|0) - |1)). 
Note that this would already complete the Fourier trans- 
form if we had only one qubit. (ii) Next, we apply to 
the first and second qubits the conditional phase shift 
\a) 2 \b)i -> e i ™ b / 2 \a) 2 \b)i (a, b = 0, 1), given explicitly by 



21 



(\ 

10 

1 

\0 e"/ 2 



B(tt/2), 



(20) 



where B(<p) is the unitary transformation 



B((j)) = exp[—icj) 



'0 





,0 -1. 



(21) 



Then we "split" the second qubit by applying to it the 
transformation A. This accounts for the contribution of 
the second most significant bit of the input x. (iii) Sim- 
ilarly, in steps i = 3,4, ... ,9 we apply the conditional 
phase shift \a)i\b)j -> e 17rab / r ~ 3 \a) l \b) :j (a,b = 0, 1), that 
is, 



B; 



(\ 

10 

1 

\0 e™l T 



B(ir/T-i), (22) 



to each pair of qubits j = 1,2, — 1, and then 

apply the transformation Aj = A to the ith qubit. Note 
that after the ith. step the first i qubits store the Fourier 
transform of the i most significant bits of x. Hence, after 
the 9th step the Fourier transform is completed: 



F 2 a — (A 9 B 9S ■ ■ ■ B. 



91 i 



(A 3 B 32 B 31 )(A 2 B 21 )(A 1 ). (23) 



Performing these operations implies also application of 
state exchanges whenever one needs to transfer the states 
of atoms i and j to a common unit cell for processing. A 
list of control commands (SCk sequences) corresponding 
to Eqs. (p~8|) and (p3|) can be written straightforwardly. 



VI. CONCLUSION 

We have shown that quantum devices with exponen- 
tially large Hilbert space dimension can be efficiently 
controlled, provided they are assembled from completely 
controllable unit cells in an architecture that is optimized 
for the specific function they should perform. The unit 
cell can be constructed from simple quantum objects of 
arbitrary physical nature: two-level atoms, nuclear spins, 
rotating molecules, quantum dots, etc. This allows to 
optimize critical properties such as coherence time and 
control precision for practical realizations. The only re- 
quirement is that the unit cell could be put under non- 
holonomic control, i.e., that it could be sufficiently per- 
turbed to have unconstrained dynamics. This ensures 
that the cell can be fully controlled and made perform 
any desired operation. 

As a concrete example, we have considered a quantum 
system of 2™ levels, composed of n two-level atoms that 
are coupled by dipole-dipole interactions. The atoms are 
grouped into unit cells, each consisting of three nearly 
resonant atoms. Each cell is controlled with two time- 
dependent perturbations: a static electric field and an 
electromagnetic field nearly resonant with the atoms. 
We have shown that any unitary transformation in the 
2 3 = 8 dimensional Hilbert space of the cell can be ef- 
fected within a few control cycles, each comprising 64 
applications of the perturbations with values fixed ac- 
cording to a non-holonomic control scheme. In particu- 
lar, the Toffoli-gate transformation on the cell regarded 
as a 3-qubit register and any permutation of the three 
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qubits can be performed. We have given two examples 
of function-specific devices that can be assembled from 
such cells: (i) By arranging the cells in a ternary tree- 
like structure, we obtain a device that can perform ef- 
ficient quantum computations on n qubits: any unitary 
transformation on any three qubits can be effected within 
order of log 3 n control cycles. We have described a toy 
device that can perform computations on 9 qubits, in- 
cluding, for example, the discrete Fourier transform, (ii) 
When the atoms are arranged in a planar lattice struc- 
ture, where at each control cycle the closest neighboring 
atoms are differently grouped in triads, we can simulate 
various linear and non-linear lattice models of quantum 
fields with time dependent parameters. 
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(A5) 



In our control scheme the Toffoli-gate transformation can 
be effected on the unit cell by repeating 8 times the trans- 
formation U e /s = exp(— iirli.Tof f /$>) (e = tt), which is 
directly attainable: U(t = 64T) = U e/8 (see Fig. |(d)). 



The authors are grateful to M. Gromov for discussions 
and stimulating remarks. G.H. acknowledges support 
from the Foundation for Fundamental Research on Mat- 
ter (FOM), which is financially supported by the Nether- 
lands Organization for Scientific Research (NWO). 



APPENDIX A: 

The Toffoli-gate transformation is the unitary trans- 
formation on three qubits, 

U T off\x 2 )\xi)\xo) = |a? 2 )|a!i)|aJoXOR(xi ANDasa)), 

(Al) 

which corresponds to the three-bit classical logic gate, 



■T-2 
Xl 
Xo 



X2 
Xi 



x XOR (xi AND x 2 ), 



(A2) 



introduced by Toffoli as a universal gate for classical 
reversible computation \22\. It acts as a permutation 
of the computational basis states, \x) = \x2)\xi)\xq) , 
x = J2r=o Xr ^ T ~ 0, 1, . . . , 7, given by the unitary ma- 
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This matrix can be presented as 

Uroff = exp(~iTTH T off), 
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